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IJSV/ EEICTIOKAL RESISTANCE LAW JOE SMOOTH PLATES* 
By P. Schultz-Crunow 



SUMMART 



Prom measurements in the free "boundary layer of a 
plate the laws governing the velocity distribution and a 
now resistance law are derived which, by increasing 
Reynolds number Rej^ afford lower resistance values than 
the logarithmic law* The transverse velocities, the shear- 
ing stress, and the mixing path profiles were also defined. 

IIJTRODUCTIOK 



The application of the logarithmic laws of velocity 
distribution for turbulent pipe flow (reference l) to the 
free friction layer had afforded a resistance law which, 
after minor changes of the experimental constants contained 
in the velocity distribution laws, could be brought into 
satisfactory agreement with the plate drag measurements, up 
to the highest Re^ numbers. 

nevertheless, this application constitutes no more 
than an approximation for the still xinknown velocity dis- 
tribution in the free friction layer, for there is no cogent 
necessity for an identical velocity distribution in the 
pipe and on the plate; one may only surmise that they dif- 
fer slightly from one another. Aside from that the plate 
drag measurements (references 2 and 3) are not completely 
satisfactory, since they were achieved with comparatively . 
small test plates on which the assumption of plane flow is 
net only in the neighborhood of the plate leading edge or 
else obtained on not completely hydraul ical ly smooth plates 
(reference 4). It therefore seemed desirable to explore 
the velocity distribution in the free friction layer and 
to check the drag measurements. 



*"]ileues Heibungswiderstandsgesetz fur glatte Flatten," 

Luf tf ahrtf orschung, vol. 17, no. 8, Aug. 20, 1940, pp. 239- 

46. 
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Ifotation 

X distance from plate leading edge - coordinate of 
plate length 

y wall distance 

z coordinate of plate width 

I mixing path 

u velocity in direction x 

u mean velocity in friction layer 

U flow velocity 

V velocity in direction y 

T shearing stress 

Tg wall shearing stress 

p density 



p rate of shearing stress 

6 friction layer thickness 
8 

5*= / ^1 - dy displacement thickness 

0 ,6 

- ■ I ^ ( 1 - dy momentum thickness 
0 'B_lsa_ local drag coefficient 

/ T dr ■ ' 

0 0 ■ 

Cj = total drag coefficient 

P ^ 

I " - 

u kinematic viscosity 
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Hex = 
Hes = 




i. Beynolds numbers 



SXPBRIMBlTTiiL SET-UP 



The measTxrements were to "be made in the tunnel so 
as to insure greater accuracy, while at the same time the 
conditions on free surfaces were to be preserved; i.e.* 
no pressure decrease, and a free friction layer, which 
means that the friction layer is followed by a parallel 
floV with constant velocity distribution. A newly erected 
experimental arrangeir;ent for air operation was enployed. 
It consisted essentially of a blower-operated tunnel of 
rectangular section (fig. l), the lower, horizontally 
placed vjall of which carried the surface to be explored. 
The aspect ratio of the cross section was chosen with a 
view to plane flow in the median zone of the horizontal 
sides. The tunnel height was so chosen that the opposite 
friction layers were always kept separate by the nuclear 
flow, that is, the zone of uniformly distributed speed U 
in figure 1, and thus produced free friction layers. The 
upper wall of the tunnel was hinged and adjustable so that 
any prescribed pressure distribution and, for our purposes, 
also a pressure equal to the outside pressure could be ob- 
tained to within 1/20 millimeter alcohol accuracy. The 
tunnel height thereby increased in flow direction accord- 
ing to the proportional growth of the displacement thick- 
ness at the walls. This made the conditions in the tunnel 
the same as on free surfaces. 

That the flow in the median zone of the test plate 
is in fact plane can be seen on the velocity profiles in 
figure 2 plotted against the logarithm of the wall resist- 
ance y, as recorded in a plane at right angles to flow 
direction and experimental wall, in the tunnel center 
(z = 0) and 250 millimeters to the left (z = -250) and 
right (z = 250) of it as viewed in flow direction. The 
intermediate profiles, which were also recorded and lie in 
the same range of scattering, have been omitted for reasons 
of clarity. 

To make sure that the transitional region was defi- 
nitely situated on the plate leading edge, this edge was 
greatly curved, and, to insure the formation of a new 
friction layer on it, a slot was provided below this edge 
through v/hich the friction layer of the' blower chamber 
wall exhausted (fig, l) , 
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RESISTANCE MEASUREMENT 



There were two ways of measuring the friotional re- 
sistance. The momentum loss ~ / u(U - •u.)dy could "be 

dx,/ 

determined from the measurement of the dynamic pressure 

over wall distance y at different distances x from 
the plate leading edge and the wall shearing stress Tq , 
then computed from the momentum equation. 

8 



u(U - u)dy (l) 




Or the friction could he weighed directly on a rectangular 
test plate 'counted movably in a sector of the principal 
plate. The latter method proved more accurate. A siirilar 
method had already "been used elsewhere (reference 3). The 
arrangement is illustrated in figure 3. The test plate 
rests on an arm in flexure pivots I rotatable ahout a 
vertical axis. The weight is taken up hy a float S. 
Aside from the moment of the friction force an opposite 
moment is applied in the hinge hy a wire D stressed in 
torsion. The wire can he twisted with the hand wheel un- 
til it balances the moment of the friction force and the 
scale arm is in the neutral setting, which can he read 
optically. The torque for the related twist was calihrat- 
ed so that the friotional force and hence the local resist- 
ance could be ascertained. This arrangement is practical 
only in our case, where the pressure in the tunnel is the 
same as in the outside space. In any other case the slots 
in the Sector necessary for the free movement of the test 
plate manifest flows which produce uneven suction and pres- 
sure on the plate edges and falsify the measurement. 

Three surfaces were explored. Surface A was a 
huilt-up plywood plate of 25 millimeters thickness; surface 
B, the puttied, polished, and lacquered surface of a high- 
pov-rer metal airplane; and surface C consisted of separate, 
1.5 millimeter thick metal panels arranged as in figure 4, 
flush riveted on a' 1.5 millimeter thick metal plate. The 
whole was fastened with countersunk wood screws to a braced 
plywood plate and, without prior puttying, given a smooth 
camouflage coat. 
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Surface A .- The moinentum, displacement, and friction- 
layer thickness obtained fron the velocity measurements 
are reproduced in figure 5, The next quest ion -.was 
whether the zero point of the diagram could he permitted 
to coincide with the plate leading edge, hecause the lead- 
ing edge of the test plate was of finite thickness and 
was rounded off. Only a minor correction was disclosed 
hy recalculation of the plate length from the friction- 
layer thickness ohtained in the foremost test section 
(x = 0.25 m) hy the old power law (reference 5), itfhich 
at small Reynolds ntushers corresponds to reality quite 
well . 



The momentum thicl-ness d can be tied to the drag 
coefficients throxigh equation (l). It affords with the 
equation for ^: 

Cf ' = 2 ^ (2a) 
^ dx. 

Cf = 2 ~ (2h) 

so that the drag coefficients can he determined from the 
measured momentum thickness. But this determination was 
found to he not accurate enough. The scale measurements 
shown in figure 6 are much more accurate. A curve I in- 
cltLded for comparison reproduces the resistance law (ref- 
erence l) obtained with the logarithmic velocity distri- 
bution law for flow in pipes, curve II the old power law 
(reference 5) and curve III an empirical system of formu- 
las (reference J5.) , 

0.242 ,„ V 

(3a) 



°f 



lOgdlCj. Cf ) 



0,558 Cj, 

Cf' = 1 — (3b) 

0.558 + ZiTo^ 

obtained from all the plate measurements made up to that 
time. 

It may be mentioned that (3a) and the so-called 
Schlichting approximate formula (reference l) 

0.455 

= ' (4) 

f / - v2.5a 

Clog Re } 
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which refers to curve I are practically -in accord within 
our range of experimentation. 

Figure 6 manifests a systematic departure of the 
equalizing curve 17 of our measurements from curve I to 
the extent that the new measurements give lower drag 
coefficients at log Re^ = 7,2 for instance, the values 
are 4 percent lower. At any rate, surface A proved hy- 
draulically smooth under the conditions of our measrue- 
ments* 

As to the accuracy of these measurements it should 
be added that the test surface was slightly wavy, since a 
technical surface was involved. Even though it did not 
show up as roughness, the plate nevertheless manifested 
areas where the edges of the test plate Stood out or back 
a little, thus inducing slot flows. But the effect of 
these flows could be nullified by the design of the lead- 
ing and trailing edges of the plate shown in figure 3. A 
0,1-milliraeter vertical displacement of the plate on the 
scale produced a 1.5 percent change in force. Since the 
plate height could be accurately adjusted to 0.1 milli- 
meter, this percentage is the degree of the test accuracy 

On surfaces B and C the accuracy was substantially 
less, First, because the test length was only 2 meters 
instead of 6 meters; then, no suitable pressure orifices 
could be provided on the surfaces without causing inter- 
ference; hence zero pressure could not be reliably ob- 
tained; and lastly, the edges of the plate could not, 
without destroying their surface, be sharpened in the 
requisite manner so as to minimize the effect of the slot 
flo\tf properly. Thus a 0. l-millimeter vertical displace- 
ment of the test plate had already yielded a 6-percent 
change in force, and so it may be stated merely that the 
test valties scattered around the power law and that the 
surfaces are likewise hydraui ically smooth within the ex- 
plored range (to Re^ = 10''); 



VELOCITY DISTRIBUTION 



The laws of velocity distribution in the free fric- 
tion layer are derived from the same points of view as 
for pipe flow (reference l) . There are two such laws 
(reference 7), one for wall proximity, the other for the 
remaining zone of the friction layer (reference l) . The 
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first resulted from the fact that in wall proximity the 
wall distance is the characteristic quantity for the 
flov/ attitude. It is expressed with 

For not too close wall proximity it is 



= A log ^ + B (5h) 

where a and B are., uni vers al constants and defined from 
the test. The second law stems from the argiiment that the 
s iin T ■> T tude of the flow at greater wall distances is solely 
defined ''oj the pipe radius ana the frluLion layer thi<?T»-np.Ra 
8, respectively. Here it affords 




On the pipe fg can he fairlj^ approximated from (5"b) ; 
hence the wall proximity law may he extrapolated over the 
pipe radius. The velocity measurements in figures 7 and 
8, in fact, disclose that the test points in wall proxiiaity 
and in v/all distance, respectively, form one curve. This 
is particularly -so In figure 8 across the entire test 
range with -exception of the test point closest to the wall. 
In this plot the test points near the wall terminate in a 
straight line, a, sign that the linear wall law (Sh) is 
applicable, as is plainly seen in figure 7. The constants 
here assume the values 

Ap = 5.93, Bp = 4.07 

subscript P denoting the plate. The pipe tests on the 
other hand (reference l) had given 

Ah = 5.74, Bh = 5.46* 



*T'he5e values correspond to average values of smaller and 
larger Heynolds numbers usad as a basis for the determina- 
tion of the logarithmic resistance law. For large Reynolds 
numbers the values Aj^ = 5.84 and = 5.52 are better. 
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In the calculation of the logarithmic law of plate 
resistance (reference l) with equation (Sh) these constants 
were modified to 

Aj' = 5.85, Bj^' .= 5,56 

so as to afford agreement with Kempf *8 experiments (refer- 
ence 3). Therefore, discounting the minor discrepancy 
in the experimental constants the same wall proximity lav/ 
as in the pipe is applicable. However, figure 8 shows in 
the wall distance even more than in pipe flow a distinct 
departure from this law, to the extent that the test points 
deflect downward from the straight line in order to reach 
the atscissa hefore the straight line. In this instance 
the application of the pipe flow to the plate affords 
therefore a less good approximation, for the straight line 

U - u 

gives ■ ' ■ — =0 at a too large ahscissa value and, since 

V* 

the friction layer thickness is defined with this value, 
this thickness was computed a little too high, the momen- 
tum loss too great and hence a drag a little too high. 
This explains the discrepancy of our drag measurements 
from the logarithmic resistance law. 

Jigure 8 shows, as has "been mentioned, a universal 
relationship with y/6 nearly throughout the entire test 
range* Only the test point nearest the wall diverges and 
thereby manifests the viscosity effect in the immediate 
proximity of the wall. The range of the y/6 relation 
in the individual test series fluctuates between 0.0039 <. 

^ 0.02, to which values log — « 1.2 corresponds 

to subscript' i denoting the limit of validity. In the 
presence of the smallness of 7\ I ^ it is questionable 
whether it is permissible to continue the straight line 
in figure 8 as far as the wall. The question gains in 
importance in the subsequent determination of the resist- 
ance ilaw, because of the required integration over y. 
It is therefore checked on the integrals 
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To check these values against the test, "by transforming 
with (6) to 

Jo 

-~ was determined with the compensating curve IV of fig- 

ure 6, according to the relation ( ) ~ Sc^'. The 

integrals subsequently ottained in figure 5 are given in 
the t?>T<=*- They vary very little from the constant mean 
values which are in complete agreement with ilm jjraTicusly 
secured values, so that the extrapolation as far as the 
wall is definitely permis sitle. Vfith these data the re- 
lations for the displacement and momentum thickness follorfs 
from (?) at 

= 2.36 J 7 = T- - Cf" 

o 5 5 



RESISTANCE LAW 



At first it might Tie thought expedient without regard 
to the estahlished velocity distribution to retain the 
logarithmic resistance law and to "bring it into accord 
with the compensating cuve IV of figure 6 "by subsequent 
change of the constants A^, Bg, But that is not possible, 
as it merely would afford a parallel shift of the old 
curve I; whereas curve IV has a distinctly different in- 
clination and hence would require a rotation of curve I 
also. 

The new plate resistance law is therefore derived 
from equation (l) on the basis of our universal velocity 
distribution (fig. 8), The function fg may be intro- 
duced in (l) , since its extrapolation as far as the wall 
has proved permissible. It affords with (?) 



o o 
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and, with the above values for the integrals, after replac- 
ing X, 8:Bej^, Heg : 

(T/Tt:h(--^---(T/)] <«' 

V* 

As for the two unknown functions Reg for which 

a relation is necessary, it may be stated that the veloc- 
ity distribution can be tied in with y/6 as well as 

yv* 

— as previously indicated on equations (5) and (6). 
Specifically the relation 



-rr = a, + b In 



— -j-^ = aa - b In ^ 



(9a) 
(9b) 



affords from its identity the desired relation 



or 



1. •. v*8 
aj. + as + b log — ^ 



-(-- ^) 

Ee^ ■ ^ e^^'' ■ ^ (10) 



with a = a + ag. 

Kow the tangent of the velocity distribution in a 
poiiit of the wall adjacent viscid fluid layer is to be ex- 
pressed by (9), a and b being temporarily left unde- 
fined and these constants later so chosen that the resist- 
ance law obtained agrees in two point, that is, at two Rej. 
numbers with the experiment, that is, the compensating 
curve I? in figure 6, So. rather than make an approximate 
assumption about the velocity distribution in direct wall 
proximity, a certain tangent is picked from out of the 




respectively, the detailed 
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Equation (lO) written in (s) gives the differential 
equat ion 



iff.. 



)r 



(—) -TT^iT^' [3,34 rjj- - 21.4(^j J| 

for ^ and hence for c-.', since = -/? Abtrevi- 



U 

ated it affords with 



Sr= T,; Hex = t = P 



d i 

III 



Id 



whence hy differentiation 
1 



= T,' e° 



1 / 

~ ! 3.34 

h \ 



31.4 \ ^ 21.4 
TJ / Tl 



and ty separation of variables 



:^(3.34 - 21.4 Ti) + 21.4 



dTi 



The integration affords 

t - to^ = [3.34 - (?, X 3.34 d + 21.4)ti ^ 

+ 2 X 3.34 h^ + 21.4 h + 21.4] 



i/e integrate from = 0. Then ii^ = 0 (v* = oo), 

whence the lower integrf.tion limit on the right-hand side 
is 2 X 3.341)2 + 21.4 "b + 21.4. The value from the numeri- 
cal evaluation ranged at around 100, vrhile the upper in- 
tegration limit lies at 10'^; in consequence of which the 
lovjer limit can he suisinarily disregarded. Then the solu- 
tion readsJ 

In t = ^ - P + ln[3.34 - (6.64 h + 21.4)ti 

+ 6.64 h^ + 2l.4(b + l)] (ll) 
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which expresses our new resistance law, namely, the rela- 
tionship 'between c^ • and Re. Now the constants h and 
P are so defined that (ll) agrees with the experimental 
curve in two points* 

This occurs when 

h = 2.895; P = 1.25 

and so defines our tangent in the viscid fluid layer, with 
the exception of its position, since the velocity distri- 
hution in this immediate wall proximity is without effect 
on the momentum integral. 

With the previously cited values of b and P it 
afforded agreement not only in the two points with the 
compensating curve IV of figure 6, "but practically through- 
out the intermediate zone as well, thus furnishing an add- 
ed proof of the universal validity of the obtained veloc?? 
ity distribution laws, and so permits an extrapolation of 
our resistance curve to any Reynolds number as exemplified 
in figure 9 as far as Re^^ = 2 x 10^, The logarithmic law 
is also shown for comparison. 

The total resistance coefficient c^ was numerically 
defined from figure 9 with the relation 



X 




0 



readily derived from (s). It is reproduced in figure 10, 
along with the logarithmic law and the old power law. The 
departure of the logarithmic law is quite noticeable at 
the highest Ke^. 

The approximate formulas for the two new resistance 
laws covering the range 10 ® < Rojj < 10 are as follows* 

, 0.370 
°f " (log Re,,) "-^"^ 

g _ 0.427 

• ^ (-0.407 + log Re^) ^'^"^ 

They rcorrespond to the form of approximate formula proposed 
by Prandtl (reference l) 
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(log He- + p) 



n 



with, which three points can oe made to agree exactly and 
v^hich was used also in the derivation of (4); p itself 
was found to "be practically equal to zero in the deriva- 
tion of the Cj* formula, as in the derivation of (4). 



3ESAE STEESS AUD MIXIIIG PATH DISTRIBUTIOU 



The shear stress follows from the application of the 
m on en tun equation 



5 — / u dy + u V = ^ 



TABLE I.- EVALUATIOlT OP VELOCITY MEASUREMENT 



TAverage values: fs d ^ = 3 . 32 , ( f 2 ) ^ d = 21.4] 

" ■ no 00 



X 

(m 


) 


1 

U i 

(c E /s ) 


log 
,(He,) 


V* 

'u 


5 

(cm) 


S* 

(cm) 

t 


(cm) 


1 

/f2dy/8 


1 

/fs^dy/S 


^0 
(cm/s) 


0. 


5 


1 

1 

193? : 


5.831 


0.0444 


1. 


30 


0.187 


0.131 


3.20 


21.9 




1 




1935! 6.132 

I 


.0414 


2. 


24 


.306 


.222 


3.30 


21.9 


3.96 


1. 


5 


1942 


6. 30 8 


.0398 


3. 


07 


.410 


.301 


. 3. 36 


22. 3 


3.67 


2. 


5 


1941 


6, 531 


.0381 


4. 


6 2 


.634 


.446 


3.30 


20.6 


S.20 


3. 


2 


1941 


. 538 


.0373 


5. 


62 


.700 


.538 


3. 34 


20,8 


3.30 


3. 


9 


1938 


5.723 


.0367 


6. 


60 


.810 


. 624 


3. 35 


21.0 


3.45 


5. 


3 


1940 


6.857 


.0358 




40 


1.022 


.796 


3.40 


21.0 
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Blf f erentiation leaves 





0 

The continuity equation 



V = - 

also exists. For the determination of — the velocity 

ax 

distribution was plotted in figures 11 and 13 against dif- 
ferent wall distances. The curve of is shown in fig- 

ures 13, and 14, and that of the transverse velocity in 
figure 15; Vq denotes the raaximum transverse velocity 
which ocqurs at the edge of the friction layer. The sepa- 
rate Vg values, compiled in the table disclose a maximum 
Vq of about 2 percent of the flow velocity. Further 
evaluation gives the shear stress profiles of figure 16, 
concerning which it is stated that the equation of motion 
of a fluid particle past the plate (no pressure gradient) 
reads 

u~ 4.V — - - — ~ 
9x Sy p ay 

But, at the wall, the fluid adheres whence ^ =0 at that 

point. This condition for the initial tangent is not al- 
ways satisfactory complied with on the reproduced shear 
stress profiles because of inaccuracies in the evaluation 
which are unavoidable on account of the differentiation of 
the v(x) curve. At the most, the wall shear stresses ob- 
tained from this momentum interpretation are 5 percent 
greater than those obtained by balance measurement accord- 
ing to which the reliability of the evaluation can be gaged. 

The mixing length is obtained from 
and shown in 'figure 17. It now approaches at - = 0,4, 
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the constant mean value = 0.072. The relation for the 

5 

wall tangent, 1 = 0,43y, does not vary apprecialDly from 
the relation = 0,43y* ascertained on the pipe. A 

difference must exist "by reason of the different shear 
stress profiles. 

Translation by J. Vaaier, 
national Advisory Committee 
for Aeronautics. 
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Pigs. 1,2.3 



Section /IS 




Figure 1.- Experimental setup. 
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Figure 2,- Velocity profiles on the 

experimental wall; T in mm. 



Figure 3.- Arrangement of 
test plate and 
lialance. 

H, band wheel 
D, wire 

T, spring hinge 
S. float 



KACA Technical Memorandum No. 986 



Figs. 4,5 




Figure 4,- Surface C; flush riveted; arraoigement of rivets and screws, 

a, flush rivets 

1), countersxmk wood screws, not covered. 




Figure 5,- Momenttun displacement and friction layer 
thickneas 5,6*, ^ on surface A (imilt-up 
Ic—iinated plate) , U = 19.4ai/8, 
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6.- Local drag coefficient of stirface A, 

I, calculated fromKenipf' 8 pipe test (logarithmic 

resistance lav) , 
II, old power fonnula Cf '0,0576 He -1/5 
III, £clioenherr'« empirical fonmila, 
IV, conipenBatixig curve through test points. 
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Figure 9.- New resistance law for smooth 
plate, Cf' curve. 



I. l^e?/ resistance law. 
II. Losarithinic resistance law. 
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Figare 10.- New resistance law for smooth 
plate, cf cui've. 
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Hgs. 13, 14, 15. 




Figure 15,- Transverse velocity 
U = 19.4 m/s. 
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